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Abstract 

In this work, we discuss the neutrino masses and mixings as the realization of an 
S3 flavour permutational symmetry in two models, namely the Standard Model and an 
extension of the Standard Model with three Higgs doublets. In the S3 Standard Model, 
mass matrices of the same generic form are obtained for the neutrinos and charged leptons 
when the S3 flavour symmetry is broken sequentially according to the chain S31, <8> S3R D 
gdiag j n minimal S3-symmetric extension of the Standard Model, the S3 

symmetry is left unbroken, and the concept of flavour is extended to the Higgs sector by 
introducing in the theory three Higgs fields which are SU(2) doublets. In both models, 
the mass matrices of the neutrino and charged leptons are reparametrized in terms of their 
eigenvalues, and exact, explicit analytical expressions for the neutrino mixing angles as 
functions of the masses of neutrinos and charged leptons are obtained. In the case of the S3 
Standard Model, from a x 2 fit of the theoretical expressions of the lepton mixing matrix to 
the values extracted from experiment, the numerical values of the neutrino mixing angles 
are obtained in excellent agreement with experimental data. In the S3 extension of the 
Standard Model, if two of the right handed neutrinos masses are degenerate, the reactor 
and atmospheric mixing angles are determined by the masses of the charged leptons, 
yielding #23 in excellent agreement with experimental data, and #13 different from zero but 
very small. If the masses of the three right handed neutrinos are assumed to be different, 
then it is possible to get #13 also in very good agreement with experimental data. We 
also show the branching ratios of some selected flavour changing neutral currents (FCNC) 
process as well as the contribution of the exchange of a neutral flavour changing scalar to 
the anomaly of the magnetic moment of the muon. 
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1 Introduction 



The observation of flavour oscillations of solar, atmospheric, reactor, and accelerator neutrinos 
established that they have non-vanishing masses and mix among themselves, much like the 
quarks do p~M20] . In these observations and experiments, the differences of the squared neutrino 
masses as well as the neutrino mixing angles are measured. These discoveries brought out very 
forcibly the need of extending the Standard Model (SM) in order to accommodate in the 
theory the new data on neutrino physics in a consistent way, that would allow for a unified and 
systematic treatment of the observed hierarchy of masses and mixings of all fermions. At the 
same time, the number of free parameters in the extended form of the SM had to be drastically 
reduced in order to give predictive power to the theory. These two seemingly contradictory 
demands are met by a flavour symmetry under which the families transform in a non-trivial 
fashion. The observed pattern of neutrino mixing and, in particular, the non vanishing and 
sizable value of the reactor mixing angle strongly suggests a flavour permutational symmetry 
S3. 

The result of a combined analysis of all available neutrino oscillation data, including the 
recent results from long-baseline — > v e searches at the Tokai to Kamioka (T2K) [21] and 
Double CHOOZ experiments [22] as well as the Main Injector Neutrino Oscillation Search 
(MINOS) experiment [23], give the following values for the differences of the squared neutrino 
masses and the mixing angles in the lepton mixing matrix, U PMNS , at la confidence level 
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the upper (lower) row corresponds to inverted (normal) neutrino mass hierarchy, see also Gon- 
zalez Garcia et al [25l[26] and J. F. W. Valle et al [2?]. In fact, from the three angles needed to 
describe the mixing of the neutrinos, the least known is 9 l3 . A global analysis of the T2K [21], 
MINOS [23] and CHOOZ [22] experiments yielded a non- vanishing value for the reactor mixing 
angle j27J|28]. Recently, the Daya Bay [29] and RENO experiments [30] found the following 
values for the reactor neutrino mixing angle: sin 2 29[ 3 = 0.092±0.016 (stat) ±0.005 (syst) which 
is equivalent to 9[ 3 ~ 8.8°±0.8° at 5.2 a level, and sin 2 2^ 3 = 0.113±0.013 (stat) ±0.019 (syst) 
which is equivalent to 9[ 3 ~ 9.8° at 4.9 a level [30] . 

In the last ten years, important theoretical advances have been made in the understanding 
of the mechanisms for fermion mass generation and flavour mixing. The imposition of flavour 



0.52 ±0.06 [ 0.016 H _ . 



2 



symmetries in the Standard Model and its extensions strongly constrains the number of free 
parameters in the Yukawa couplings and gives rise to special forms of the fermions mass matrices 
with few free parameters and a number of texture zeroes [3TH34"] . For a recent review of flavour 
symmetry models see [351136] . 

In the case of the Minimal S^-Invariant Extension of the Standard Model [37H13] . the 
concept of flavour and generations is extended to the Higgs sector in such a way that all the 
matter fields - Higgs, quarks, and lepton fields, including the right-handed neutrino fields- 
have three species and transform under the flavour symmetry group as the three dimensional 
representation 1 © 2 of the permutational group S3. A model with more than one Higgs 
SU (2) doublet has tree level flavour changing neutral currents whose exchange may give rise to 
lepton flavour violating processes and may also contribute to the anomalous magnetic moment 
of the muon. An effective test of the phenomenological success of the model is obtained by 
verifying that all flavour changing neutral current processes and the magnetic anomaly of the 
muon, computed in the ^-Invariant extended form of the Standard Model, agree with the 
experimental values. 

Another important application of the permutational group S3 is the classification of physi- 
cally equivalent mass matrices. Different mass matrices with texture zeroes located in different 
positions have exactly the same physical content if they are related by a similarity transforma- 
tion in flavour space [3H[32j|34l|44]. ^ ^ e i nvariari ts of the mass matrix are to be preserved, 
the elements on the diagonal can only exchange positions in the diagonal while the off diagonal 
elements can only exchange positions off the diagonal. So, the transformations matrices that 
define the similarity classes are the six elements of the three dimensional real representation of 
the group of permutations S 3 [33||4"5|. 

The left-handed Majorana neutrinos naturally acquire their small masses through the type 
I seesaw mechanism 

M^M^M^M^,, (3) 

where M„ D and Mj, s denote the Dirac and right handed Majorana neutrino mass matrices, 
respectively. 

The mass matrices are diagonalized by bi-unitary transformations as 

U^MiUjij ; = diag(m i i,m i2 ,m i3 ), and Uj'M^U,, = diag(m w , m V2 , m V3 ), (4) 

where i = d,u, e. The entries in the diagonal matrices may be complex, so the physical masses 
are their absolute values. 

The lepton flavour mixing matrices V PMJVS arise from the mismatch between diagonalization 
of the mass matrices of charged leptons and left-handed neutrinos, 

V PMNS = U^K, (5) 
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where K is the diagonal matrix of the Majorana phase factors. Therefore, in order to obtain 
the unitary matrices appearing in eq. (0) and get predictions for the flavour mixing angles and 
CP violating phases, we should specify the mass matrices. Also, in the case of three neutrino 
mixing there are three CP violation rephasing invariants jlEj, associated with the three CP 
violating phases present in the V PMJVS matrix. The rephasing invariant related to the Dirac 
phase, analogous to the Jarlskog invariant in the quark sector, is given by: 

Ji = %rn [KTOKa^i] ■ (6) 

The rephasing invariant J\ controls the magnitude of CP violation effects in neutrino oscillations 
and is a directly observable quantity. The other two rephasing invariants associated with 
the two Majorana phases in the V PMJVS matrix, can be chosen as: S\ = 3m [K1K3] anc ^ 
S2 = Sm [V^V^j]. These rephasing invariants are not uniquely defined, but J/, Si and S2 are 
relevant for the definition of the effective Majorana neutrino mass, m ee , in the neutrinoless 
double beta decay. 

In the standard PDG parametrization jU], the entries in the lepton mixing matrix are 
parametrized in terms of the mixing angles and phases. Thus, the mixing angles are related to 
the observable moduli of lepton V PMJVS through the relations: 

sin p 12 — Yq^f ) sm y 23 — l-iVeaf' y i3 — l^l • {<) 

The neutrino oscillations do not provide information about either the absolute mass scale 
or their nature, this is, if neutrinos are Dirac or Majorana particles [15] . Thus, one of the most 
fundamental problems in neutrinos physics is the question of the nature of massive neutrinos. 
A direct way to reveal the nature of massive neutrinos is to investigate processes in which 
the total lepton number is not conserved [19]. The matrix elements for these processes are 
proportional to the effective Majorana neutrino masses, whose magnitudes squared are 

\(mi)\ 2 = Y?j=x rn l j \ v ij\ 4 + 2 Y?j<k m v 3 m »k \ v ij\ 2 \Vik\ 2 cos2 (wij - wit) , I — e, ji, r, (8) 



where m Vj are the neutrino Majorana masses, Vij are the elements of the lepton mixing ma- 
trix, wij = arg{V/j}; this term includes phases of both types, Dirac and Majorana. The 

half-life of the neutrinoless double beta decay (0^2/3 decay) can be expressed as 

G Qv |M ^| 2 I (m ee ) | 2 , where G 0v is a phase factor, \Mq v \ is the isotope specific nuclear decay matrix 
element and | (m ee ) \ is the magnitude of the effective Majorana neutrino masses defined in eq.flH]). 
Current experiment data yield only an upper bound for this quantity |(m ee )| < 0.3 eV [50|l5T]. 
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2 An S3 flavour symmetry in the SM 



In analogy with the work of A. Mondragon and E. Rodriguez Jauregui on the S3 flavour 
symmetry in the quark sector of the Standard Model [Hl[52], we will start by assuming the 
obvious, the one Higgs boson in the SM is an SU(2)l doublet and, since it has no flavour, it 
can only be accommodated in a singlet representation of S3. The mass term in the Lagrangian, 
obtained by taking the vacuum expectation value of the Higgs field in the quark and lepton 
Yukawa couplings, gives rise to mass matrices M d , M u , M; and [33] 153]; 

£y = qd,LM d q djR + q u ,LM u q UtR + L L MiL R + U L M V (u L ) c + h.c. (9) 

If it is assumed that S3 is an exact symmetry of the model, these mass matrices give mass 
only to the one fermion in each family that is assigned to the singlet representation of S3. 
Hence, in a symmetry adapted basis, all entries in these matrices should vanish except for the 
one at the third row and third column. There is no mismatch between the diagonalization 
of the mass matrices of the charged leptons and neutrinos, or d- and u- type quarks, and, 
consequently there is no mixing of the flavour indices. Therefore, we propose, along with many 
other authors [211 EH E3 SUES], that the texture zeroes of the mass matrices of quarks and 
leptons are the result of a flavour permutational symmetry S 3 and its spontaneous or explicit 
breaking. In particular, the permutational S3 flavour symmetry and its sequential explicit 
breaking allows us to justify using the same generic form for the mass matrices of all Dirac 
fermions [551I53], this form is conventionally called a two texture zeroes form. Some reasons 
to propose the validity of a matrix with two texture zeroes as a universal form for the mass 
matrices of all Dirac fermions in the theory are the following: 

1. The idea of S3 flavour symmetry and its explicit breaking has been succesfully realized 
as a mass matrix with two texture zeroes in the quark sector to describe the strong mass 
hierarchy of up and down type quarks [331 15l3] - |56] . Also, the numerical values of the 
mixing matrices of the quarks determined in this framework are in good agreement with 
the experimental data [531133] . 

2. Since the mass spectrum of the charged leptons exhibits a hierarchy similar to the quark's 
one, it would be natural to consider the same S3 symmetry and its explicit breaking to 
justify the use of the same generic form with two texture zeroes for the charged lepton 
mass matrix. 

3. As for the Dirac neutrinos, we have no direct information about the absolute values or the 
relative values of the Dirac neutrino masses, but the mass matrix with two texture zeroes 
can be obtained from an SO (10) Grand Unified Theory which describes well the data on 
masses and mixings of Majorana neutrinos [57H59] . Furthermore, from supersymmetry 
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arguments, it would be sensible to assume that the Dirac neutrinos have a mass hierarchy 
similar to that of the u-quarks and it would be natural to take for the Dirac neutrino 
mass matrix also a matrix with two texture zeroes. 



2.1 Mass matrices from the breaking of S^l (g) S^r 

Some authors have pointed out that realistic Dirac fermion mass matrices result from the flavour 
permutational symmetry S^l ® S&R an d its spontaneous or explicit breaking according to the 
chain: S 3L x S 3R D S 3 diag D S 2L x S 2R D S 2 diag [3I1E31I311EZIIS1I521I531ISD]. The group S 3 treats 
three objects symmetrically, while the structure 1 © 2 of its 3 x 3 matrix representations treats 
the generations differently and adapts itself readily to the hierarchical nature of the fermion 
mass spectra. As explained above, under exact S 3 l ® Szr symmetry, the mass spectrum for 
either quark sector (up or down quarks) or leptonic sector (charged leptons or Dirac neutrinos) 
consists of one massive particle in a singlet irreducible representation and a pair of massless 
particles in a doublet irreducible representation of S 3 l <8> S 3R . In order to be more precise, in 
the case of exact S 3 l <8> S 3R symmetry, and assuming that there is only one Higgs boson in 
the theory, this SU(2)l doublet can only be in a singlet representation of S3, that is, it is a 
scalar with respect to the S3 transformations. Hence, the corresponding mass matrices, M i3 , 
are invariant with respect to a permutation of the family (columns) and flavour (rows ) indices, 
and take the form 

M£° = ^r^l > i = u,d,l,u D 




the subindex W stands for weak basis. In order to make explicit the assignment of particles 
to irreducible representations of 5*3, it will be convenient to make a change of basis from the 
weak basis to a symmetry adapted or hiererchical basis by means of the unitary matrix that 
diagonalizes the matrix Mj3, 



where 



C^mTu (11) 



U = -== I -V3 1 y/2 I and mJ J = m i3 I (12) 




V3 


1 V2\ 




-V3 


1 y/2 


and 





-2 y/2 J 







In the Standard Model with the S3 symmetry, masses for the first and second families are 
generated if we add the terms and Ma to M^. The term Mj2 breaks the permutational 
symmetry S 3 l <8> ^3^ down to S 2 l ®> S 2 r and mixes the singlet and doublet representation of S3, 
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while the term Mji transform as the mixed symmetry term of the doublet complex tensorial 
representation of the S^ ms diagonal subgroup of S^l <%> S^r. Thus, 

An 1A2 -Ai - iA i2 \ 
-iAa ai ^ (13) 

ft ft -2ft ) w 

Finally, adding the mass matrices eqs. f lT0|) - (|T3|) . we get the mass matrix Mj in the weak basis. 
Then, in a symmetry adapted basis Mj takes the form 

/ Ai \ 

M, = m i3 J A* Bi Ci J , i = u, d, I, v D (14) 
V Ci DiJ H 

where Ai = |Aj|e l< ^, Bi = — Aj + Si and Di — 1 — Si. From the strong hierarchy of the masses 
of the Dirac fermions, m i3 >> m i2 > mn, we expect 1 — Si to be very close to unity. The 
Hermitian mass matrix (TT41 may be written in terms of a real symmetric matrix Mj and a 
diagonal matrix of phases Pj = diag [l, e 1 ^, e z ^~\ as follows: 

Mi = PlMiP u i = u,d,l,v D . (15) 

Each possible symmetry breaking pattern is now characterized by the flavour symmetry 
breaking parameter Z^ 2 , which is defined as the ratio ZAI 2 = |^'| 23 . This ratio measures the 
mixing of the singlet and doublet irreducible representations of 5*3. The small parameter Si is 
a function of the flavour symmetry breaking parameter Z^ 2 |33ll34ll4T| . 

Thus, we obtain a universal form for the mass matrices of all Dirac fermions in the theory. 
But in the Standard Model and its extensions considering a mass term for left-handed neutrinos 
purely of Dirac nature is not theoretically favored, because it cannot explain naturally why 
neutrinos are much lighter than the charged leptons. Thus, we assume that the neutrinos have 
Majorana masses and acquire their small masses through the type I seesaw mechanism . 

2.2 Classification of mass matrices with texture zeroes in equiva- 
lence classes 

In this section we make a classification of mass matrices with texture zeroes in terms of similarity 
clases [33]- The similarity classes are defined as follows: Two matrices M and M' are similar if 
there exists an invertible matrix T such that M' = TMT -1 or M' = T _1 MT. The equivalence 
classes associated with a similarity transformation are called similarity classes. Another way to 
see the similarity classes is that the matrices that satisfy the similarity transformation have the 
same invariants: trace, determinant and X — \ (Tr {M 2 } — Tr {M} 2 ) . Therefore, all matrices in 
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a class of similarity have the same eigenvalues, since all have the same characteristic polynomial, 
given by Af - Tr {M 2 } Af - x\ ~ det {M} = 0. 

Now, from the most general form of the symmetric and Hermitian mass matrices of 3 x 3: 





M = \ a b c and M — la* b c , (16) 



we can see that only six of the nine elements of these matrices are independent of each other. 
Therefore, the similarity transformation is realized as the permutation of the six independent 
elements in the nine entries of the mass matrices. But if we want to preserve the invariants, the 
elements on the diagonal can only exchange positions on the diagonal, while the off-diagonal 
elements can only exchange positions outside the diagonal. Thus, all these operations reduce 
to the permutations of three objects. So it is natural to propose as transformation matrices T 
in the similarity clases, the six elements of the real representation of the group of permutations 
5*3 which are: 



T(A ) = 1 \ ,T(A 1 )=\ 1 , T(A 2 ) 




T (A 3 ) = 1 , T (A 4 ) = 1 , T (As) 





(17) 



Then, we get the classification of mass matrices with texture zeroes, which is shown in the 
table 1. In this table, the "*" and "x" denote an arbitrary non- vanishing matrix element on 
the diagonal and off-diagonal entries, respectively. We recall the rule for counting the texture 
zeroes in a mass matrix: two texture zeroes off-diagonal counts as one zero, while one on the 
diagonal counts as one [31] . 



2.3 The mass matrix for left-handed neutrinos 

The left-handed Majorana neutrinos acquire their small masses through the type I seesaw 
mechanism , eq. fl3]). The form of M„ D is given in eq. (TH|) . which is a Hermitian matrix with 
two texture zeroes of class I. From the conjecture of a universal S3 flavour symmetry in a unified 
treatment of all fermions, it is natural to take for also a matrix with two texture zeroes 
of class I, non Hermitian but symmetric. Let us further assume that the phases in the entries 
of M UR may be factorized out as [33l|3l]: = RM^R, where R = diag [e"^ c , e^ c , l] with 
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Table 1: Matrix with two (left) and one (right) texture zeroes 




arg { c u } and 







a. 



■r 







M Vr = I a VR \b UR \ \c Vr \ 
\c Ur \ d VR 

Then, the mass matrix of the left-handed Majorana neutrinos has also the same generic 
with two texture zeroes of class I: 



form 



M, 








d Vj 



(19) 



where 



a, R ' a »L ~ d VR ' C " L ~ d VR + \ C " D e d » R J 



(20) 



The elements Oj, and d v are real, while b v and are complex. Therefore, the form of 
mass matrices with two texture zeroes is invariant under the action of the seesaw mechanism 
of type I ^HE21EaiSI]. 

It may also be noticed that, if we set b VR = or/and c Vr = 0, the resulting expression for 
M„ still has two texture zeroes. Therefore, M,, also has two texture zeroes when M,, has 

L ' L U R 
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two, three or four texture zeroes (the last two cases are called Fritzsch textures) [33j[3l]. In 
fact, the information on the total number of texture zeroes present in the mass matrix M v , 
can be obtained from the elements (2, 2) and (2, 3) 23 of the left-handed neutrinos mass matrix. 

From the previous analysis, the matrix M„ has two non-ignorable phases which are 0i = 
argjo^} and (p 2 = ax s{°i/ L }- However, to describe the phenomenology of neutrinos masses 
and mixing, only one phase in is required. Therefore, without loss of generality, we may 
choose 0i = 202 = 2y? [33j[34] . In this case the analysis simplifies, since the phases in M Vi may 
be factorized out as 

M^ = QM^Q, (21) 

where Q is a diagonal matrix of phases Q = diag [e~ lip , e lip , l] and M„ is a real symetric matrix. 
Then, the matrix M v , can be diagonalized by a unitary matrix of the form U„ = QOj,, where 
O v is an real orthogonal matrix that diagonalizes the real symetric matrix M„ . 



2.4 Mass matrix as function of the fermion masses 



The real symetric matrix Mj, eqs. (fT5l) and (l2Tj) . may be brought to diagonal form by the 
transformation, Mj = Ojdiag {mji, m i2 , mj 3 } Of, where the mis are the eigenvalues of Mj and 
Oj is a real orthogonal matrix, with i = u, d, I, v L . Now, computing the invariants of the real 



symetric matrix Mj, we may express the real parameters aj, bi, Ci and d{ occuring in eqs. (TT5 
and (I2ip in terms of the mass eigenvalues [331IM] as: 

2 m iX m i2 m i?> 



mi + m i2 + mz - di 



(di - run) (di - m i2 ) (m i3 - di) 



(22) 



di di 

From the condition that aj, bi, q and di are real, we determine the allowed region of di [62J. In 
other words, all elements of the matrix Mj must be real and depending on which eigenvalue is 
chosen as negative, the parameter di must satisfy one of the following conditions: 

Normal hierarchy Inverted hierarchy 





, > di> m i2 , 


for 


mji = - 


m 


mi 


, > di > mji, 


for 


m2 = - 


m 


m. 


> di> ma, 


for 


mz = - 


m 



m i2 > di> m i3 , for mji 
mji > di> m i3 , for m i2 
m i2 > di> m iu for m i3 



\mi\, 
\mt2l 
\m3\- 



In this way, for a normal hierarchy and taking m i2 = — \rrii 2 \, we get the Mj matrix 
u, d, I, u L ), reparametrized in terms of its eigenvalues and the parameter Si as 



/ 







1-Si 







1-Si 



mi - m2 + fa 



filf-, 



i2 



1 - Si 



(23) 
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where fh a = m i2 = fa = l- fha - S { and f i2 = 1 + m i2 - S { . 

The small parameters Si, which are also functions of the mass ratios and the flavour symme- 
try breaking parameters Z X J 2 , are obtained as the solution of the cubic equation (1 — Si)(fha — 
ma + Si) 2 Zi = Sifnfi2- The solution of the cubic equation that vanishes when Zi vanishes may 
be written as 

A - Zi ("^2 ~ fhg) 2 

l ~ Zi + 1 W t {Zi) ' [ZV 

where Wi (Z) is the product of the two roots of the cubic equation that do not vanish when Zj 
vanishes 



with 



and 



Wi(Z) 



pf + 2qf + 2qJpf + q, 



Qi + \/pf + Q 2 



+ 



\Pi\ + 
i \ 



Pi + Qi 



p? + 2^-2 ft ^p? + , 
[Zi (2 {m i2 - mji) + I' 



+ 



Pi 



+ {m i2 ~ mil) + 2) + | (Zj (2 (m i2 - ma) + 1) + [m i2 - fha) + 2) 2 , 



1 Zi(Zi(2(mj2-mq)+l)+mj2-"til+2) 2 , [Zj (m i2 -rag ) (mjg -rag +2) (l+m^) (1-mq )] 



1 [Zi(mi2-ma)(mi2-mq+2)(l+mi2)(l-m.ii)](Zi(2(rni2-rn.ii)+l)+rni2-m,ii+2) 
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(Zi+iy 



(25) 



1 (Zi(2(mi2-mii)+l)+mj2-raii+2) 

27 (Xi+l) a ' 

The allowed values for the parameters Si are in the following range < Si < 1 — fha. 

Now, the entries in the real orthogonal matrix O that diagonalize the matrix Mj, may be 
expressed as 







i 




1 




1 


ma fn 


2 


fhiifia 
T> i2 


2 


T>i3 


2 






1 




1 




1 



Oi 



ma(l-5i)fn 


2 


mi2{l-Si)fi 2 


2 


'0-Si)Si 


2 




1 


V i2 


1 


T> i3 


1 



V 





2 


mizfiiSi 


2 


filfi2 






T> i2 







(26) 



/ 



where, T> a = (1 - ^)(ma + m i2 )(l - ma), V i2 = (1 - <y(m»i + m i2 )(l + m i2 ), and V i3 
(1 - - fha)(l + m i2 ). 



2.5 Mixing Matrices as Functions of the Fermion Masses 

The unitary matrices \J u ,i occurring in the definition of V PMJVS , eq. (jSJ), may be written in polar 
form as = P^/O^i. In this expresion, P^/ is the diagonal matrix of phases appearing in 
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the two texutre zeroes mass matrix (JT51) . Then, the lepton mixing matrix takes the form 

V th = OfP^O.K, (27) 

PMNS 1 v ' v > 

where p(^~') = diag [l, e i$1 , e** 2 ] is the diagonal matrix of the Dirac phases, with $ x = 2tp — <pi 
and $2 = <p — 4>i- The real orthogonal matrices Ut i are defined in eq. (}2~6l) . Thus, the mixing 
matrix V PM whose entries are explicit function of the masses of the leptons has the following 
form [M] : 

Kmns = V* V%e^ , (28) 

where 

v£ = + (v / F¥ r »^ + V^NUe^) , 

= ~V IZvJ: 1 + Vl£ (V(l-*)(l-W«/,ie^ + v/WW^) , 

= + (Va-^Cl-Wza/^ + VWn/.ie^) , (29) 

\r th _ / mem^m^m^S^ / 5i6„(l-6i)(l-6„) jfy / fnfnfvifa 

in these expresions the f^s, and 2Vs are 

/i/(J)l = (1 - «V(e) - > = (1 - £/(!)) (™*i(e) + "VMX 1 _ "Vi(e))> 

/i/(i)2 = (l + fh^fjt) - 8 v (j,)) , V v{ i) 2 = (1 - (m^fe) + mn,o*))(l + ^Kj(m))> ( 30 ) 

A/(I)3 = (1 - ^(0)( 1 ~ + "Mm))> 

where fh Vl t e ) = m " l{e) ; and m V2 u) = l^ 2 -^!. In the quark sector, the elements of the mixing 

matrix V CKM may also be expressed as functions of the quark masses ratios, the resulting 
expressions are similar to the expressions obtained above for the elements of the matrix V PMJVS , 
for more details see 



12 



2.6 The x 2 fit for the Lepton Mixing Matrix 

We made a \ 2 fit °f the theoretical expressions for the modulii of the entries of the lepton 
mixing matrix \{V PMNS )ij\ given in eq. ( |29i) to the values extracted from experiment as given 
by Gonzalez- Garcia and Maltoni [14]. The computation was made using the following values 
for the charged lepton masses [47] : 

m e = 0.51099 MeV, m„ = 105.6583 MeV, m T = 1776.82 MeV. (31) 

We took for the masses of the left-handed Majorana neutrinos a normal hierarchy. This allows 
us to write the left-handed Majorana neutrino mass ratios in terms of the neutrino squared 
mass differences and the neutrino mass m U3 in the following form: 

_ n = ^ = ^7^. (32) 

The numerical values of the neutrino squared mass differences were obtained from the global 
analysis of the experimental data on neutrino oscillations given in [H], and we left the mass 
free parameter of the \ 2 fit- Also, the parameters 8 e , 8 V , $i and $2 were left as free 
parameters to be varied. Hence, in this \ 2 fit we have four degrees of freedom. 
From the best values obtained for m Ua and the experimental values of Am| 2 and Am^, we 
obtained the following best values for the neutrino masses at la: 

m vi = (3.22±g;|J) x 10" 3 eV, m V2 = (9.10±^;f|) x 10" 3 eV, m U3 = (4.92±° f 2 ) x 10~ 2 eV. (33) 

The resulting best values of the symmetry breaking parameters are 61 — (6 ± 2.98) x 10~ 2 and 
5 V = 0.522to' ( j ) 2; an d the best values of the Dirac CP violating phases are $1 = (270 ± 15)° 
and $2 = (180 ± 10)°. The best values obtained for the modulii of the entries of the PMNS 
mixing matrix are given in the following expression 

PMNS 

The value of the rephasing invariant related to the Dirac phase is j\ = (1.8 ± 0.6) x 10~ 2 . In the 
absence of experimental information about the Majorana phases Pi and /?2, the two rephasing 
invariants Si and S2 associated with the two Majorana phases in the V PMJVS matrix, cannot 
be determined from experimental values. Therefore, in order to make a numerical estimate 
of the Majorana phases, we maximized the rephasing invariants Si and S%, thus obtaining a 
numerical value for the Majorana phases fli and fa- Then, the maximum values of the rephasing 
invariants are S^ ax = -4.9 x 10~ 2 and S% ax = 3.4 x 10" 2 , with ft = -1.4° and f3 2 = 77°. In this 



0.8204; 
0.3748; 
0.4345- 



-0.008 
-0.010 
-0.018 
-0.031 
-0.024 
-0.020 



0.5616; 
0.6280; 
0.5388 



-0.012 

-0.014 
-0.019 
-0.010 
-0.022 
-0.024 



0.118llffi 
0.6819 ±0.025 



0.7216 



-0.024 
-0.027 



(34) 
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numerical analysis, the minimum value of the x 2 ■, corresponding to the best fit, is x 2 = 0.288 
and the resulting value of x 2 f° r degree of freedom is ^j^f = 0.075. All numerical results of the 
fit are in very good agreement with the values of the moduli of the entries in the matrix ~V PMNS 
as given in ref. [T4"] . 



2.7 The Mixing Angles 

The theoretical expressions for the neutrino mixing angles as functions of the charged lepton 
and neutrino mass ratios are obtained from eqs. ([7]), when the theoretical expressions for the 
modulii of the entries in the PMNS mixing matrix are substituted for |Vy| in the right hand 
side of eqs. ([7]). If we keep only terms at leading order, we obtain: 



sin 2 C « {m l /T Z - \ J , (35) 



. 2 nl th ^ <*v + Sefv2 ~ VS u S e fu2 COS ($ 1 - $^ 
sin (7 2 3 ~ / _ \ 













m M J 


(1- 


H m U2 ) 



(36) 



sin 2 9$ « V^, il. (37) 

( 1 +i?)( 1 +^ 2 



"2 , 



From eqs. ( 130]) we have that = 1 + — The expressions quoted above are written in 
terms of the ratios of the lepton masses, defined in eq. (1231) . When the well known values of 
the charged lepton masses, the values of the neutrino masses, eq. ( 155|) . the values of the delta 
parameters and the Dirac CP violating phases obtained from the x 2 fit in the lepton sector, 
are inserted in eqs. (|35]) -( 13~T|) . we obtain the following numerical values for the mixing angles 



C = (34.43t°;|)° , C = (43.601®° , C = (6.80lg^)° , (38) 

which are in very good agreement with the latest experimental data pT^l I2TH2T] . We may 
conclude that: 

1. The strong mass hierarchy of the Dirac fermions produces small and very small mass 
ratios of charged leptons. Then, the Dirac fermions mass hierarchy is reflected in a 
similar hierarchy of small and very small Dirac fermion flavour mixing angles [34] . 



2. The normal seesaw mechanism type I which gives very small masses to the left-handed 
Majorana neutrinos with relatively large values of the neutrino mass ratio m Ul /m U2 allows 
for large 9[ 2 and 9 l 23 mixing angles (see eqs. fl35|) -f l37|) ) . 
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In the quark sector, the mixing angles may also be expressed as functions of the quark masses 
ratios, the resulting expressions are similar to the expressions obtained above for the lepton 
mixing angles, for more details see [31]. 

2.8 The effective Major ana masses 

The theoretical expression for the squared magnitude of the effective Majorana neutrino mass 
of the electron neutrino is: 

\(m ee ) th \ 2 = m 2 Vl \V e i\ A + m 2 U2 \V e2 \ A + m 2 3 |V e3 | 4 + 2m vl m V2 \V e i\ 2 |V e 2| 2 cos (w el - w e2 ) ^ 
+2m Ul m U3 | V e i| 2 \V e3 \ 2 cos (w el - w e3 ) + 2m U2 m U3 \V e2 \ 2 \V e3 \ 2 cos (w e2 - w e3 ) , 

where w e i = argjV^j}. In a similar way, the theoretical expression for the squared magnitude 
of the effective Majorana neutrino mass of the muon neutrino is: 

\{m^) th \ 2 = m 2 ui | V e i| 4 + ml 2 |Ve 2 | 4 + m 2 3 \ V e3 \ A + 2m Ul m U2 \ V e i\ 2 |Ve 2 | 2 cos (w el - w e2 ) 

+2m Ul m v . 3 | Vei| 2 |V e3 | 2 cos (w el - w e3 ) + 2m U2 m v . 3 \V e2 \ 2 \V e3 \ 2 cos (w e2 - w e3 ) , 

where u> M , = arglV^,}. Expressions for the squared magnitude of effective Majorana masses 
I (jn ee ) th I and I (m^) th I with only terms of leading order are given in Barranco et al [31] . From 



these expressions and the numerical values of the neutrinos masses given in eq. f[33|) . we obtain 
the following expressions for the effective Majorana masses with the phases as free parameters: 

|(m ee ) th | 2 fa {9.41 + 8.29 cos(l° - 2ft) + 4.3cos(l° - 2w e3 ) + 4.31 cos 2(ft - w e3 )} x 1(T 6 eV 2 , 

(41) 

where w e3 = arctan { ° f 5 lZo%~Sl }• Similarly, 

\{m^) th \ 2 w {4.8 + 0.17 cos 2(44° -w^) + 1.8cos2(^ 2 - w^)} x 10~ 4 eV 2 , (42) 

where w, 2 « arctan { HKflff" } and ^ - arctan { jfgfe^ }. 

In order to make a numerical estimate of the effective Majorana neutrinos masses |(m ee )| 
and |(m wt )|, we used the following values for the Majorana phases ft = —1.4° and ft = 77° 
obtained by maximizing the rephasing invariants S\ and S 2 . Then, the numerical value of the 
effective Majorana neutrino masses are: | (m ee ) th \ fa 4.6 x 10~ 3 eV, | (rn^) 1 ^ fa 2.1 x 10~ 2 eV. 
The numerical value of | (m^) | obtained in this theoretical scheme is at least two orders of 
magnitude smaller than the current experimental upper bound. Therefore, it will be very 
difficult to determine with good precision the numerical value of the effective mass of Majorana 
neutrinos, since, at present, this value could be of the same order magnitude as the error bars 
of the nuclear matrix elements. 
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3 A Minimal invariant Extension of the Standard Model 
with three Higgs bosons 

In the Standard Model analogous fermions in different generations have identical couplings to 
all gauge bosons of the strong, weak and electromagnetic interactions. Prior to electroweak 
symmetry breaking, the Lagrangian is chiral and invariant with respect to permutations of 
the left and right fermionic fields, that is, it is invariant under the action of the group of 
permutations acting on the flavour indices of the matter fields. 

Since the Standard Model has only one Higgs SU(2)l doublet, which can only be an S3 
singlet, it can only give mass to the quark or charged leptons in the S 3 singlet representation, 
one in each family, without breaking the S3 symmetry. Hence, in order to impose S3 as a 
fundamental symmetry, unbroken at the Fermi scale, we are led to extend the Higgs sector of 
the theory [37]. The quark, lepton and Higgs fields are Q T = (ul, di), Ur, d R , L T = (u L , e^), e R , 
ur, and H, in an obvious notation. All these fields have three species, and we assume that each 
one forms a reducible representation Is © 2 of the S3 group. The doublets carry capital indices 
/ and J, which run from 1 to 2 and the singlets are denoted by Q3, U3R, d 3 R, L 3 , e3R, V3R and 
Hs- Note that the subscript 3 denotes the singlet representation and not the third generation. 
The most general renormalizable Yukawa interactions for the lepton sector of this model are 
given by Ly = Cy E + £y„ , where 



£y b = -YfL I H s e I R-Y. 3 e L3Hse3 R -Y^[L I K I jH 1 e JR + L I r] I jH2ejR} 
-YiL 3 H I e I R-Y 5 e L I H I e3R+ h.c, 

C Yv = -Y»L^ia 2 )H* s v I R-Y»L3±ia 2 )H* s v 3R -Y»L 3 (ia 2 )H* I v I R 

-Y 2 U [ L I K IJ (ia 2 )HlvjR + L I 7] I j{ia 2 )H^vjR ] - Y£ L I (ia 2 )H* I V3 R + h.c, 



and 

' 1 \ , / 1 



(43) 
(44) 

\ J and i={o -1 )■ ^ 

Furthermore, we add to the Lagrangian the Majorana mass terms for the right-handed neutri- 
nos, 

£ M = -^CM VR u R7 (46) 

where C is the charge conjugation matrix and = diag{Mi, M 2 , M3} is the mass matrix 
for the right-handed neutrinos. In the following we will consider two cases: i) when the first 
two masses are degenerate, M\ = M 2 7^ M3, ii) when all masses are different. 

Due to the presence of three Higgs fields, the Higgs potential of the S3 invariant extension 
of the Standard Model is more complicated than that of the Standard Model [38l[63l|6l] . This 
potential was first analyzed by Pakvasa and Sugawara [51] who found that in addition to the S3 
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Table 2: Z 2 assignment in the leptonic sector. 





+ 


Hs, v 3R 


Hi, L 3 , Li, e^R, em, vm 



symmetry, it has an accidental permutational symmetry S' 2 : Hi •f-)- H2, which is not a subgroup 
of the flavour group S3. In this communication, we will assume that the vacuum respects the 
accidental S' 2 symmetry of the Higgs potential and that (Hi) = (H 2 ). With these assumptions, 
the Yukawa interactions, eqs. (I43p - (j44p yield mass matrices, for all fermions in the theory, of 
the general form [37] 

+ A*2 A*2 A*5 \ 

y,i 111 - fi 2 /i 5 • (47) 

In principle, all entries in the mass matrices can be complex since there is no restriction 
coming from the flavour symmetry S3. Therefore, there are 4 x 5 = 20 complex parameters 
in the mass matrices, which should be compared with 4 x 9 = 36 of the SM, including the 
Majorana masses of the left-handed neutrinos. 



3.1 The mass matrices in the leptonic sector and Z2 symmetry 

A further reduction of the number of parameters in the leptonic sector may be achieved by means 
of an Abelian Z 2 symmetry [37]. A possible set of charge assignments of Z 2 , compatible with 
the experimental data on masses and mixings in the leptonic sector is given in Table [2j These 
Z 2 assignments forbid the following Yukawa couplings Y± = Y% = = Y^ = 0. Therefore, 
the corresponding entries in the mass matrices vanish, i.e., \i\ — /i| = and \i\ = ^ = 0. 
The mass matrix of the charged leptons takes the form 

(fa h h \ 
h -fa fa ■ (48) 
fa fa / 

The resulting expression for M e , reparametrized in terms of its eigenvalues and written to order 
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m M m e /m 2 ) and x A = (m e /m^) 

( 



is 



M e ps m T 



V2 s/T+x 2 



V2VTT2 



fh e (l+x 2 ) „jS e 



V2VT- 



s/2 Vl+x 2 



m e (l+x 2 ) „jS e 



1 

V2 



V2 



1+x 2 — m? 



\ 







(49) 



y y/l+x 2 — rh 2 ^ y/l+x 2 — m 2 J 

This approximation is numerically exact up to order 10 -9 in units of the r mass. Notice that 
this matrix has no free parameters other than the Dirac phase S e . 

The unitary matrix U e L that enters in the definition of the mixing matrix, V PMJVS , is 
calculated from Ug L M e MjU e £ = diag (m^m^mj), where m e , and m T are the masses of 
the charged leptons, written to the same order of magnitude as M e is XJ e i = P e O e , where 
P e = diag {l, 1, e iSe } and 



/ -4 



{l+2m 2 ^+Ax 2 +fhj l +2m 2 ,) 
x/2"" ^l+m 2 +5x 2 -fnf L -mf l +m, 2 +12x 4 



{l-2fh 2 l +fh'l-2m 2 ) 



{1+Ax 2 



-2ml) 



-5x 2 



hj.-ml+fh 2 +l2x i 



\f2 ^l-4m 2 +x 2 +6mf l -4rhf l -5m 2 

1 (l-a^+m^) 
V2 ^l-4m 2 +x 2 +6mf l -4rhf l -5m 2 



1 

V2 



^l+2x 2 —rh'^—rh 2 (1+m^+x 2 — 2r~n 2 ) 
\ yjl+m 2 ^+5x 2 -fhf l -m^+fh 2 +12x i 



| {l+x 2 -™ 2 -2m 2 )^/l+2x 2 -m 2 -m 2 y r T+^rh e m^ 
^l—Am^+x 2 +6mj l —4m^—5rh 2 ^J\ J ^x 2 —m l L 



in this expression, as before, m M = m^/m T: m e = m e /m T and x = m e /m^ 



(50) 



3.1.1 Neutrino masses and mixings with two degenerate masses of the right- 
handed neutrinos 

According to the Z 2 selection rule the mass matrix of the Dirac neutrinos takes the form 

[331321130! 

/ & & \ 

M UD = \ & -f4 , (51) 

and considering the following form to the mass matrix of right-handed neutrinos = 
diag {Mi, Mi, M3}, the mass matrix for the left-handed Majorana neutrinos, M„ L , obtained 
from the seesaw mechanism type I is, 

2p\pl 

M Vi = I 2(p») 2 I , (52) 

1{P\? + (P3T 




18 



where p v 2 = (/i^/M^ 2 , p\ = (jifi/M?* and pi = (pQ/M^- M x and M 3 are the masses of the 
right handed neutrinos. The complex symmetric neutrino mass matrix M Vi may be brought to 
a diagonal form by the transformation U^M^U,, = diag (Im^Je^ 1 , \m U2 \e l< ^ 2 , \m U3 \e % ^ , where 
XJ U is the matrix that diagonalizes the matrix Mj, Mj, L . This allows us to reparametrize the 
matrices and XJ U in terms of the complex neutrino masses, which have the following form: 



,1/2 



rV2. 



■in 



M 







m 



\/(rn V3 - m Ul ){m V2 - m„ 3 )e t5v 




(53) 



s/ (m U3 - m Vl )(m V2 - m U3 )e l5u 
and U„ = P^O^ where = diag {l, 1, e tSv } and 

/ 

/ cos?7 sm?7 u \ 



[m 



m 



1'2 



m V3 )e 



cos rj sin t] 
1 
— sin r) cos r) 



V 





/ m "3 







\ 
















1 




/ rra„ 3 —m Ul 




-m„ 3 





/ 


m„ 2 — m„j 


Y "i„ 2 





(54) 



The unitarity of Uj, constrains sin 77 to be real and thus | sin 77 1 < 1, this condition fixes the 
phases <p\ and (f) 2 as \rn Ul \ sin0i = \m U2 \ sin 02 = l^i^l sin^. The only free parameters in the 
matrices and XJ U are the phase <j) u , implicit in m Ul , wn V2 and m U3 , and the Dirac phase 8 U . 

The neutrino mixing matrix V PMJVS , is the product U^U^K, where K = diag (l, e ia , e 4/3 ) 
is the diagonal matrix of the Majorana phase factors. Therefore, the theoretical mixing matrix 



V 



PMNS ' 



is given by 



/ 0\ x cos rj + 31 sin rje lSl 0\ x sin 77 — 0\ x cos rje zSl — Ofi \ 



V 



-Of 2 cos r\ + 0\ 2 sin r]t %bl —0\ 2 sin 77 — 0\ 2 cos rje l51 O 
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K, 



(55) 



y Of 3 cos T) — O33 sin r\e %&l Of 3 sin 77 + O33 cos Tye* 5 ' O23 / 
where cos 77 and sin 77 are given in eq. f l54"|) . Of^ are given in eq. (l5"0]h and 5/ = 5^ — <5 e . 
The mixing angles 

The magnitudes of the reactor and atmospheric mixing angles, 9 1 ^ and 9\\ , are determined by 
the masses of the charged leptons only. Keeping only terms of order (ml/mfj and (m At /m T ) 4 , 
we get [20] 

1 _ (l+4x 2 -m4) 



smt% « —a; 



'23 



/)/* 

sin 6>? 



1 !+^ 2 -2m2 +7 n4 
\/2 v /l-4m2+a; 2 +6m4 



(56) 
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The magnitude of the solar angle depends on the charged lepton and neutrino masses as well 
as the Dirac and Major ana phases 



/ 



tan0 



12 



1 -2^cos«J, 

O31 1 



V 



2^cos<5, 
O31 



m„ 3 


- m Vl 


m V2 


- m Vz 




- m vz 


m U3 


- m Ul 



On 
On 



m 



+ 



m 



m 



(Qn) 

\o 3 i J m V3 - m Vl ) 



(57) 



The dependence of tan#f 2 011 the Dirac phase 5i, see ( )57|) . is very weak, since O31 ~ 1 but 
On ~ l/"\/2(^e/' m /x)- Hence, we may neglect it when comparing (|57|) with the data on neutrino 
mixings. The dependence of tan Q 1 ^ on the phase 4> u and the physical masses of the neutrinos 
enters through the ratio of the neutrino mass differences, it can be made explicit with the help 
of the unitarity constraint on XJ U as 



m 



(\m V2 


2 _ 


m„ 3 


2 SiB 2 ^) 1 / 2 - 


m V3 \ 


COS 4>v 




2 _ 


m V3 


2 sin 2 0,) 1 /2 + 


m U3 \ 


COS (j> v 



Similarly, the Majorana phases are given by sin 2a = sin(0i — <p 2 ) and sin 2/3 = sin( 
where 



sin 2a 
sin 2(3 



|m„ 3 | si\i(f> v 
\m vl \\m V2 \ 



m„„ r - m,. 



I 2 sin 2 (j) u + 



m,. 



f^(|m i , 3 | v / l-sin 2 ^ + ^ 



(58) 
(59) 



m,„ r — m,. 



I 2 sin 2 > 



A more complete and detailed discussion of the neutrino mixing matrix Vpmns an d the Majo- 
rana phases obtained in this model is given in refs. [3?J|39j|l0j|l3] . 



In the present ^-invariant extension of the Standard Model, the experimental restriction 
|Am 2 2 | < |Am 2 3 | implies an inverted neutrino mass spectrum \m V3 \ < \m Vl \ < \m U2 \ [571 HP] . 
In this model, the reactor and atmospheric mixing angles, 6[ 3 and 9 23 , are determined by the 
masses of the charged leptons only as 

sin 2 0jj = 1.1 x 10" 5 , sin 2 9% = 0.5. (60) 

Even in this simplified analysis it is clear that the S3 symmetry gives a non-vanishing reactor 
mixing angle, within the bounds of MINOS [23], albeit small. The atmospheric angle is in very 
good agreement with the recent experimental data. As can be seen from equations ( |5"7|) and 
(1581 , the solar angle is sensitive to the differences of the squared neutrino masses and the phase 
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4> v but is only weakly sensitive to the charged lepton masses. If the small terms proportional 
to On an d 0\ x are neglected in ([57]) . we obtain 



tan tf 12 = — — . (blj 



(Amf 2 + Am 2 3 + 


m„ 3 


2 cos 2 0,) 1 / 2 - 


\m V3 


COS 4>u 


(Am 2 3 + 


m 2 , , cos 


2 0,) 1 / 2 + 


m V3 \ 


cos <; 





From this equation, we may readily derive expressions for the neutrino masses in terms of 
tan Q[ 2 , cos (f> u and the differences of the squared neutrino masses 



Am 2 3 1 - tan 4 0f 2 + r 2 



2 tan C cos fa J 1 + tan 2 ^ A + tan 2 ^ + r 2 



(62) 



Im^J = ^/l 771 ^! 2 + Am 2 3 , and \m U2 \ = a/|vt1j/ 3 | 2 + Am 2 3 (l + r 2 ), where r 2 = Am 2 2 /Am 2 3 « 
3 x 10~ 2 . As r 2 << 1, the sum of the neutrino masses is 

y | m I ~ V^jjS A + 2 L ~ 2 tan 2 ^ cos + tan 4 g i£ _ tan 2 _ ( 63) 

1 2 cos 4> u tan 6^ 2 V / 

The most restrictive cosmological upper bound for this sum is \m v \ < 0.l7eV [18]. This upper 



bound and the experimentally determined values of tan 612 and Am|, give a lower bound for 
cos^j, > 0.55 or < <p u < 57°. The neutrino masses \m Vi \ assume their minimal values when 
cos0„ = 1. When cos(p u takes values in the range 0.55 < cos0 < 1, the neutrino masses change 
very slowly with cos0j,. In the absence of experimental information we will assume that cf) v 
vanishes. Hence, setting <j) u = in our formula, we find 

m Vl = 0.052 eV m U2 = 0.053 eV m Vz = 0.019 eV. (64) 

The computed sum of the neutrino masses is (X)i=i \Tn Vi \) = 0.13 eV, which is consistent with 
the cosmological upper bound [18], as expected, since we used the cosmological bound to fix 
the bound on cos(p u . The effective Majorana mass in neutrinoless double beta decay (m ee ), is 
defined in eq. (jSj). The most stringent bound on (m ee ), obtained from the analysis of the data 
collected by the Heidelberg-Moscow experiment on the process neutrinoless double beta decay 
in enriched Ge [SO], is (m ee ) < 0.3 eV. In this framework for a preliminary analysis we are 
assuming that the Majorana phases vanish, thus we get 

(m ee ) th = 0.053 eV (65) 

well below the experimental upper bound. 
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Deviation of the mixing matrix Vp MJVS from the tri-bimaximal form 

The previous results on neutrino masses and mixings weakly depend on the Dirac phase 5, 
for simplicity we will assume in this work that S = 7r/2. We may write the mixing matrix as 
follows, V 



th 



V + AV where the tri-bimaximal form V ESI is 

PMNS PMNS PMNS ' ' 



\rtri 
v PAINS 



and AV% NS 



AV e + 5t 



(V2 + 5t 



12, 



12" 



g(5t 



AV,. 



(66) 



12 J 



where g(Sti 2 ) = 1 + §<5ti2(\/2 + 5t\ 2 ). All entries in AV e are proportional to (me/m^) 2 except 
(AV e ) 13 that is proportional to (m e /m M ). 

The value for <5t 12 , which is a small parameter, fixes the scale and the origin of the neutrino 
mass matrix. If we take for 5t\ 2 the experimental central value 5ti 2 ~ —0.04, we obtain 
\m U2 \ « 0.056 eV, ImJ « 0.055 eV, and \m m \ « 0.022 eV gOl When we take for 5t 12 the 



tri-bimaximal value 5t\ 2 = 0, the neutrino masses are m 



0.0521 eV, m U2 



0.0528 eV, and 



0.0178 eV. For a detailed discussion on this subject see [H 



In both cases the S3 invariant extension of the SM predicts an inverted hierarchy. Since 
the tri-bimaximal value of 5t 12 differs from the experimental central value by less than 6% of 
tan #i2, the difference in the corresponding numerical values of the neutrino masses are not 
significative within the present experimental uncertainties. 



3.1.2 Neutrino masses and mixings when the masses of the right-handed neutrinos 
are non-degenerate 

In the minimal S'3-invariant extension of Standard Model, the Yukawa interactions and the 
S3 x Z 2 flavour symmetry yield a mass matrix for the Dirac neutrinos of the form flBTl) . The 
masses of the left-handed Majorana neutrinos, M^, are generated by the seesaw mechanism 
type I, eq. ([3]), where is the mass matrix of the right-handed neutrinos, which we take 
to be real and diagonal but non-degenerate = diag(Mi, M 2 , M3). Then, the mass matrix 
lSA UL takes the form 



M 



2(m!0 2 


2A(^) 2 




M 


M 


M 


2A( M £) 2 


2(^) 2 




M 


M 


M 






2(^) 2 | (^) 2 


M 


M 


M 1 M 3 



A 



\A h+M 2 J 



andM 



> M1M2 
'M2+M1 



(67) 
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When the first two right-handed neutrino masses are equal, the parameter A vanishes and we 
recover the expresion for given in Kubo et al [37], eq. (152"]) in the present paper, which 
leads to the results presented in the previous section [39H43] . 

Since we assumed the right-handed neutrino mass matrix to be real, the complex 

symmetric neutrino mass matrix M^ L has only three independent phase factors that come 
from the parameters /i2, /i3 and /14. Here, to simplify the analysis we will consider the case 



when argj/i^} = arg-f/i^} or 2axg{n u 4 } = arg < 2 



The general case, with three 



M 1 M 3 

independent phase factors, will be considered in detail elsewhere. 

In the case considered here the diagonalization of may be reduced to the diagonalization 
of a mass matrix with two texture zeroes discussed in section 12.31 The phase factors may be 



factored out of as 



M 



where = e^ 2 diag {l, 1, e i& ») with 5 U 



arg {^4} - ar g{/4}; 



(6* 




/'•» 



2|A* 



i/|2 
2 



^f(l-|A|), and M 



M 



A 
ABC 
C 2d 



(69) 



with A 



v Af 



2|A||^ 


i2 




M 





1^8 I 
M 3 



1 17 1 ^ 



|A|), C = (1 + |A|) and 



. As mentioned before, the diagonalization of is reduced to the diagonalization 



of the real symmetric matrix M, which is a matrix with two texture zeroes of class I [33]. Hence 
the matrix is diagonalized by a unitary matrix 



\J v = CfU„O v 



N[I] 



(70) 



In the literature, these similarity transformations are also known as weak basis transformations, 
since they leave invariant the gauge currents [45J. 

As in the case of the charged leptons, the matrices and \J U can be reparametrized in 
terms of the neutrino masses. For this we use the information that we already have about the 
diagonalization of a matrix with two texture zeroes of class I [331 1341 1441 [52] . Then, the mass 
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matrix for a normal [inverted] hierarchy in the mass spectrum takes the form 
/ Mo + rf d -L( C NlI] +A NlI] ) \ 



M 



N[I] 



d 



fi +d 



V2 



v ^= ( C N[I] + A Nll] ) \(C N[1] -A N[ 



V2 



m 



+ m U2 + m V3 - 2 (p + d) , 



(71) 



with C N[I] = y v "" ru -j an d A 

The values allowed for the parameters /i and 2d + fi are in the following ranges: m U2[1] > n > 

m„ 1[sl and m U3[2] > 2d + fi > m^,. The orthogonal matrix O*' 7 ' reparametrized in terms of the 
neutrino masses is given by 



A / ( M +M ~m vi ) (2d+n -m vm ) (m„ 3[2] - 
V 2d 



2d) 



V 2d 



^ /[-1](^3-Mq)("^2-Mq)/i / ( m -3[l] "MO ) (^0 - m -l[3] ) /2 / [~ 1] {Z ) (m„ 3 - Mq ) jf W ^ 



V 



jvTTT 



p. 



V 



[-l]2d( ( u -m„ 1 )/i 



2d(m„ 2 -n Q )f. 
^T7] 



N[I] 



(m^ 2 -M )/i/ : 



-l]2d(m„ 3 -M )/ 3 

(m„ 3 -/J )/i/. 
^Pl 



W[J] 



N[I] 



where 



N[I] 



V 



N[I] 



2d(m V2 - m Ul ) (m„ 3[1] - m Vl[3] J , 
2d(m 



2d(m U2 - m Ul ) (m Us[2] - m U2[3] J 



"3[1] m ^l[3] J I m ^3[2] m ^2[3] 



/i = (2d + n -m 



V\) 1 



(72) 



JV[7] 



N[I] 



■1] (m^ - - 2d) . 



-1] (2d + fi -m V2 ) , 

The superscripts TV and / denote the normal and inverted hierarchies respectively. 
The neutrino mixing matrix 

The neutrino mixing matrix V PMJVS , is the product U^U^K, where K is the diagonal matrix 
of the Majorana phase factors, defined by K = diag(l, e ia , e J/3 ). Now, we obtain the theoretical 
expression of the elements for the lepton mixing matrix, V pM , which is: 

v el v e2 e v e3 e 

th th ia th i/3 

v ill v /j,2 t v fi3 e 



V 



(73) 
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where 

V el - ^Z U 11 ~ U 2l e ) V e2 - mT U 12 ~ U 22 e > »Vl ~ U 31 ' 

Ks = frO™-0™e«>, C = -°iT-ft0 2 T^, ^ = O a T, (74) 

with 5i = 5 U — S e , the elements N[I] are given in the previous section. 
The Reactor Mixing Angle 

The theoretical expression for the lepton mixing angles as functions of the lepton mass ratios are 
readily obtained when the theoretical expressions for the modulii of the entries in the PMNS 
mixing matrix, given in eqs. ( 174"|) . are substituted for \Vij\ in the right hand side of eqs. (|7|). In 
a first, preliminary analysis for the reactor mixing angle 9 13 and for an inverted neutrino mass 
hierarchy (m V2 > m Vl > m U3 ) we obtain: 

■2ni „ (»o + 2 d- m vz ) (fi - m vz ) 

sm v 13 r — — . yio) 

(m Ul - m U3 ) (m U2 - m V3 ) 

Now, with the following values for the neutrino masses m U2 = 0.056 eV, m Vl = 0.053 eV and 
m U3 = 0.048 eV, and the parameter values 5i = ir/2, fi = 0.049 eV and d — 8 x 10 -5 eV, 
we get sin 2 B[ 3 « 0.029 — > 6[ 3 « 9.8°, in good agreement with experimental data [29l[30]. A 
more complete analysis, from a x 2 fit of the exact theoretical expressions for the modulii of the 
entries of the lepton mixing matrix of the |(V )^-| to the experimental values (for example 
the values given in the analysis by Gonzalez-Garcia [H]) will be considered in detail elsewhere. 



4 Flavour Changing Neutral Currents (FCNC) and g-2 

Models with more than one Higgs 577(2) doublet have tree level flavour changing neutral 
currents. In the Minimal ^-invariant Extension of the Standard Model considered here, there is 
one Higgs 577(2) doublet per generation coupling to all fermions. The flavour changing Yukawa 
couplings may be written in a flavour labelled, symmetry adapted weak basis as |^0 | I¥T | B3] 



£FCNC = ( EaL Y^ s E bR + U aL Y^ b s U bR + D aL Y° s D bR ) H° s + (E aL Y^E bR + U aL Y^U bR 

(76) 

+D aL Y° 1 D bR ) Hf + {E aL Y£ 2 E bR + U aL Yg 2 U bR + D aL Y° 2 D bR ) H% + h.c. 

The Yukawa couplings of immediate physical interest in the computation of the flavour changing 
neutral currents are those defined in the mass basis, according to Y^ 1 = Ul L Y^ !I U e R, where U e L 
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and U e R are the matrices that diagonalize the charged lepton mass matrix defined in eqs. (J4]). 
We obtain 00] 



Y 



El 



where m, 





( 2m e — \m e 






Vl 




1 

2 


and 




1 ~ 2 1 ~ 




m 


= 5.94xl0" 2 , m e = 2.< 


376 x 10 4 and x 



< ! 2 



/ 



V 



-m e 



m, 



1 ~ 2 



\x\ 



J 



(77) 

2.876 x 10 4 and x = me/m^ = 4.84 x 10 3 . All the non-diagonal 
elements are responsible for tree-level FCNC processes. If the S' 2 symmetry in the Higgs sector 
is preserved [51163], (if?) = (Hi) = v. 

The amplitude of the flavour violating process /i — > 3e, is proportional to YEY* ( 
the leptonic branching ratio, 



Br(fi -> 3e) = ff^ — 3e \ and r(/i -> 3e) 



^i,2yi,2y 

\ /ie ee / 



Mfr 



-> ei/z/) * ^ ' " 3 x 2 10 tt 3 

which is the dominant term, and the well known expression for — > euu) [67], give 



Then, 



(n 



Br{n -+ 3e) » 2(2 + tan 2 0)' 



mi 



.1/ 



H 



(79) 



where M# is the neutral Higgs involved in the process, whose mass we take as M# « 120 GeV^, 
and tan/3 = 1. We obtain Br(fi — >■ 3e) = 2.53 x 10 -16 , well below the experimental upper 
bound for this process, which is 1 x 10 _ 12 [68J. 

Similar computations give the numerical estimates of the branching ratios for some others 
flavour violating processes in the leptonic sector. These results, and the corresponding experi- 
mental upper bounds are shown in Table [3j In all cases considered, the theoretical estimations 
made in the framework of the minimal S^-invariant extension of the SM are well below the 
experimental upper bounds [4"0"ll4"Tll4"3"] . 



4.1 Muon anomalous magnetic moment 

In the minimal ^-invariant extension of the Standard Model we are considering here, the Z 2 
symmetry decouples the charged leptons from the Higgs boson in the S3 singlet representation. 
Therefore, at leading order of perturbation theory there are two neutral scalars and two neu- 
tral pseudoscalars whose exchange will contribute to the anomalous magnetic moment of the 
muon [40]. Since the heavier generations have larger flavour-changing couplings, the largest 
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Table 3: Leptonic FCNC processes, calculated with a ~ 120 GeV. 



FCNC processes 


Theoretical BR 


Experimental 
upper bound BR 


References 


r — y 


8.43 x 10" 14 


2 x 10~ 7 


B. Aubert et. al. (69] 


r — y [ie + e~ 


3.15 x 10~ 17 


2.7 x 10" 7 


B. Aubert et. al. [69J 


r — y \sq 


9.24 x 10" 15 


6.8 x 10~ 8 


B. Aubert et. al. [70] 


t — y cy 


5.22 x 10" ie 


1.1 x 10" 11 


B. Aubert et. al. [7T] 


\x —y 3e 


2.53 x 10~ 16 


1 x 10" 12 


U. Bellgardt et al. [68\ 


H — y 


2.42 x 10~ 20 


1.2 x lO" 11 


M. L. Brooks et al. [72] 



contribution comes from the heaviest charged leptons coupled to the lightest of the neutral 
Higgs bosons. 

After a straightforward calculation, with the help of ( 177)) . we may write 8a^ as 

Taking again Mh = 120 GeV and the upper bound for tan/3 = 14 gives an estimate of the 
largest possible contribution of the FCNC to the anomaly of the muon's magnetic moment 
Saff^ ~ 1.7 x 10 -10 . This number has to be compared with the difference between the ex- 
perimental value and the Standard Model prediction for the anomaly of the muon's magnetic 
moment Aa^ = a e * p — a^ M [73], whose numerical value is Aa M = (28.7 ± 9.1) x 10~ 10 , which 

( H) 

means Sa^ / Aa^ m 0.06. Hence, the contribution of the flavour changing neutral currents to 
the anomaly of the magnetic moment of the muon is smaller than or of the order of 6% of the 
discrepancy between the experimental value and the Standard Model prediction. 

5 Conclusions 

We have discussed the theory of the neutrino masses and mixings as the realization of an 
S3 flavour permutational symmetry in two models, the Standard Model with an S3 flavour 
symmetry and the minimal S'3-symmetric extension of the Standard Model, with three Higgs 
doublets. 

In the Standard Model the imposition of the non-abelian permutational symmetry 5*3 as a 
broken symmetry of flavour, leads to a unified treatment of masses and mixings of quarks and 
leptons in which the left-handed Majorana neutrinos acquire their masses via the type-I seesaw 
mechanism. The explicit sequential breaking of the S3 flavour group according to the chain 
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S3r<3>S 3 l D S^ ia9 D S^ 09 , is a sufficient condition to define a generic form for the mass matrices 
of all fermions in the theory. In a symmetry adapted or hierarchical basis, this generic form is 
characterized as a mass matrix with two texture zeroes of class I. All mass matrices are, then, 
reparametrized in terms of their eigenvalues j33l[3llllll|52] . After analytically diagonalizing the 
mass matrices, explicit analytical expressions for all entries in the neutrino mixing matrix are 
obtained as functions of the masses of the charged leptons and neutrinos and one CP-violating 
Dirac phase in very good agreement with all available experimental data including the recent 
measurements of the reactor angle #13 made by the T2K, Daya Bay and RENO experiments. 

In the minimal S^-invariant extension of the Standard Model, S3 is imposed as a funda- 
mental, exact symmetry in the matter sector. This assumption leads to extend the concept of 
flavour and generations to the Higgs sector. The fermion sector of the Standard Model is left 
unaltered. Hence, going to the irreducible representation of S3, the model has one SU (2) l Higgs 
doublet in the S3-singlet representation plus two SU (2) Higgs doublets in the two components 
of the S 3 -doublet representation. In this way, all the matter fields, quarks and lepton fields, 
the right-handed neutrino fields, and the Higgs fields, belong to the three dimensional l s © 2 
representation of the group S 3 . A well defined structure of the Yukawa couplings is obtained 
which permits the calculation of mass and mixing matrices as functions of the charged leptons 
and neutrino masses [38|HU]. The magnitudes of the three mixing angles are determined by 
the interplay of the flavour S3 symmetry, the charged lepton and neutrino mass hierarchies, 
and the seesaw mechanism. The solar mixing angle is almost insensitive to the value of the 
masses of the charged leptons, but its experimental value allowed us to fix the scale and origin 
of the neutrino mass spectrum. The numerical value of the atmospheric mixing angle, 6 l 23 > 
depends strongly on the masses of the charged leptons and is in very good agreement with 
the experiment. In this model the magnitude of the reactor mixing angle, 9 13 , is sensitive to 
the difference of the values of the masses of the first and second right-handed neutrinos. In 
the case where two of the neutrino masses are degenerate, #13 is different from zero but very 
small (101112] . Allowing for the masses to be non-degenerate gives a values for #13 in very good 
agreement with recent experimental data. Explicit expressions for the matrices of the Yukawa 
couplings of the leptonic sector, parameterized in terms of the leptons masses and the VEV's 
of the neutral Higgs bosons in the S3- doublet representation, can be obtained in this model. 
Taking for the neutral Higgses M Hl 2 a very conservative value (M Hl 2 « 120 GeV) , it is found 
that the numerical values of the branching ratios of the FCNC's in the leptonic sector are well 
below the corresponding experimental upper bounds by many orders of magnitude. The con- 
tribution of the flavour changing neutral currents to the anomaly of the magnetic moment of 
the muon is small (6%) but non-negligible |iO | I^T | B3]. 
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